Four-terminal resistance of an interacting quantum wire with weakly invasive contacts 
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We analyze the behavior of the four-terminal resistance, relative to the two-terminal resistance of 
an interacting quantum wire with an impurity, taking into account the invasiveness of the voltage 
probes. We consider a one-dimensional Luttinger model of spinless fermions for the wire. We treat 
the coupling to the voltage probes perturbatively, within the framework of non-equilibrium Green 
function techniques. Our investigation unveils the combined effect of impurities, electron-electron 
interactions and invasiveness of the probes on the possible occurrence of negative resistance. 
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INTRODUCTION 



Quantum transport in novel materials, is one of the 
most active areas of present research in condensed mat- 
ter physics^. The problems that arise are specially in- 
teresting in one-dimensional (ID) devices such as quan- 
tum wires and carbon nanotubes. In these cases the 
effect of electron-electron (e-e) interactions is crucial, 
leading to the so called Luttinger liquid (LL) behavior—, 
characterized by correlation functions which decay with 
interaction-dependent exponents'^ and a power law in the 
tunneling I ~ V characteristic curved. 

The actual nature of the resistance in a mesoscopic de- 
vice has been a central issue since the first milestones in 
the theory of quantum transport. Landauer proposed the 
famous setup to study quantum transport where a meso- 
scopic sample is placed between two reservoirs at different 
chemical potentials.— Then, Biittiker—, in agreement with 
experiments^, showed the fundamental relation G = uGq 
for the two-terminal conductance of a non-interacting 
quantum wire, being n the number of transverse chan- 
nels and Go = e^/h the universal conductance quantum. 
The remarkable consequence of this simple law is the fact 
that a purely non-interacting electronic system without 
any kind of inelastic scattering mechanism has a sizable 
resistance, which for a single channel device is as large 



as Gr\ 



13kQ. This resistive behavior is due to the 



coupling between the system and the reservoirs through 
which the driving voltage is applied. For this reason, 
this quantity is identified as the contact resistance of the 
ideal non-interacting setup. The mesoscopic community 
became then motivated towards the definition of an alter- 
native physical concept to describe the resistive behavior 
of the sample, free from the effects of the contact resis- 
tance. In another pioneering workS, a gedanken setup 
was proposed in order to sense the local voltage and the 
temperature. The main idea is to consider the meso- 
scopic system locally coupled to voltage probes or ther- 
mometers, represented by means of particle reservoirs. 
The latter have chemical potentials or temperatures that 
satisfy the condition of local electrochemical or thermal 



equilibrium with the mesoscopic system, which implies 
that the chemical potentials and temperatures of these 
systems are adjusted in order to get a vanishing elec- 
tronic and heat flows through the contacts to the central 
device. For the case of two voltage probes connected 
along the sample as in the sketch of Fig. [TJ the voltage 
drop corresponding to the chemical potential difference 
(/ii — /^2)/e, defines the four terminal resistance 



Rdt — 



el '■ 



(1) 



where / is the current flowing through the setup. This 
scheme to deflne the four terminal resistance was later 
implemented in the framework of scattering matrix the- 
ory for multiterminal setups^ in wires of non-interacting 
electrons with a single^^i^ and many impurities In 
Refs. [TollTTI it is clarifled that the inference of Ra from 
a calculation based in a two-terminal geometry and the 
original Landauer formulai^ may not always be correct, 
which stresses the importance of considering a genuine 
four-terminal setup to properly evaluate this quantity. 

Among other interesting features, for non-interacting 
systems, it was predicted that negative four terminal re- 
sistances are possible<^"— This is a consequence of the co- 
herent nature of the electronic propagation along a sam- 
ple where only elastic scattering processes with barriers 
or impurities can take place. These negative (longitu- 
dinal) resistances in ballistic structures were flrst mea- 
sured in the late eightiesi^. More recently, this effect 
was also observed in semiconducting structures'^. A bit 
later, the behavior of Ra was experimentally studied in 
carbon nanotubeai^. In this case, a negative value of this 
resistance was also observed within the low temperature 
regime. Let us mention that the four-terminal resistance, 
in the context of uncorrelated fermions, has been also 
extended to the case of time-dependent voltage probes, 
leading to the concept of four-terminal impedance^^. It 
is widely accepted that the Luttinger model of inter- 
acting electrons in ID is able to capture the main fea- 
tures observed in the transport experiments of carbon 
nanotubes^'i^. In particular, the power law behavior of 
the tunneling current as a function of the applied volt- 
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age and/or temperature predicted from Luttinger liquid 
theory has been experimentally observed in these sys- 
tems. Regarding the behavior of i?4t evaluated from a 
multiterminal setup in Luttinger liquids, the literature is 
restricted to Ref. [20- Previous estimates for this quantity 
were done on the basis of an interpretation of Landauer 
formula in a two terminal setup. This is due to the fact 
that quantum transport in multiterminals Luttinger liq- 
uids or models of interacting electrons is, in general, a 
rather challenging problem from the technical point of 
view. Besides Ref. [20I genuine multiterminal systems 
have been considered in Y-type geometrie o^^i^'^ , within 
linear response in the voltage and Hartree-Fock approx- 
imation of the interaction, respectively as well as in the 
study of the tunneling current of a quantum wire in the 
Fabry-Perot regimCfSi There are also some recent works 
on the effect of wires that are capacitively coupled to an 
additional reservoir . ^^'^^ 

In Ref. 113 we have considered the setup of Fig. [U 
where an infinite Luttinger wire with a single impurity, 
through which a current I flows as a response to an 
applied voltage V, is connected at two points to volt- 
age probes. Following the procedure of previous works 
for non-interacting electrons, we have considered^ii^^-l^^ 
non-invasive contacts between the wire and the voltage 
probes. We have shown that the voltage profile displays 
Friedel-like oscillations, as in the case of non-interacting 
electronsiSiiiiHj but modulated by an envelope display- 
ing a power law behavior as a function of the applied 
voltage or temperature, with an exponent depending on 
the electron-electron interaction strength. However, it is 
known that in the opposite limit of strong enough cou- 
pling between the mesoscopic device and the probes, in- 
elastic scattering events and classical resistive behavior 
take place21. Moreover, ideal non-invasive probes cannot 
be easily realized in experimental situations. For this 
reason, the aim of the present study is to go a step be- 
yond the assumption of non-invasive probes by consider- 
ing probes that, while still weakly coupled to the sample, 
introduce decoherence through inelastic scattering pro- 
cesses, as well as inter-probe interference effects. Among 
other interesting questions, our goal is to answer if fea- 



tures in the behavior of the four-terminal resistance de- 
termined by non-invasive probes, like Friedel oscillations, 
or a negative value of this quantity, are still possible when 
the coupling of the probes becomes invasive. We address 
these issues in the framework of non-equilibrium Green 
functions and a perturbative treatment in the coupling 
to the probes. 

The work is organized as follows. In section II, we 
present the model and the theoretical treatment to eval- 
uate i?4t. In Section III and IV, we present results for the 
clean wire, and the wire with an impurity, respectively. 
Finally, we present a summary and conclusions in Section 
V. Some technical details are presented in an Appendix. 





FIG. 1: (Color online) Sketch of the setup: A voltage V is im- 
posed on a Luttinger liquid, through the chemical potentials 
for the left and right movers: /Xr,; = /i ± V/2. Two voltage 
probes are connected at the positions xi,X2- The correspond- 
ing chemical potentials /ii_2 are fixed by the condition of zero 
current through the contacts. An impurity is located at Xb- 



II. THEORETICAL TREATMENT 



Model 



We consider the setup of FigH] As in Ref. HO, we use 
the following action to describe the full system: 



S — »S*^ 



wire ^ ^imp ^ Sj-es ~t~ Scont (2) 



Si; 



where Swire corresponds to an infinite-length Luttinger 
wire and reads 



J 



dx dt 



Mi 



r 1 r '1^1 

Tpi + ipl i{dt + dx) - fj-r i>r- g i^li^r + i>]i>i I 



(3) 



r 



with the first two terms corresponding to free spinless left 
and right movers respectively and g is the e-e interaction 
in the forward channel. We use units where ?i = fc^ = 1. 
We also take the Fermi velocity of the electrons vp = 
1 and the electronic charge e = 1. The two chemical 
potentials fj,i = fi — V/2 and fir — fi + V/2, for the left 
and right species, respectively, represent a voltage bias 



V between the left and right ends of the wire, which 
generates a current /. 

The effect of the impurity is described by a backward 
scattering term with strength A_b at a given position Xb- 



Simp ^ >^B J dxdtS{x - Xb) [e~ 



-2ikFX ji 



iili^i + h.C.J . (4) 
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We describe the voltage probes by Sres, corresponding to 
non-interacting electrons with two chiralities 



where 



+xl^[iidt + dy^) - fij Xr,}, 



(5) 



The term Scant represents the tunneling between the 
reservoirs and the wire, 



E 

X pTiikpx+k'-j, 



dxdyj WjS{x - Xj)6{yj - y°) 



h.c. 



(6) 



The upper and lower sign corresponds to I and r, respec- 
tively, while kp and fc^^ are the Fermi vectors of the wire 
and reservoirs, respectively. Note that here we are assum- 
ing that the voltage probes couple symmetrically to the 
left and right movers in the LL. This is a natural assump- 
tion in the absence of magnetic fields and spin-orbit in- 
teractions (Recall that we are considering a spinless LL). 
Since our main interest is to discuss the effects originated 
in the strength of the couplings, later, in Sections III and 
IV, as a simplifying hypothesis, we will consider the same 
coupling for both probes (wi = 1^2 ) , although it is known 
that an asymmetrical coupling (wi ^ W2) is sufficient to 
produce a negative resistance^i^. 

The tunneling current from the probes to the wire is 



(7) 



Gl;,fi^{x,yf,t,t') = i < xl^{yj,t')^c.{x,t) >, (8) 



is the lesser Green function involving degrees of freedom 
of the wire and reservoirs. 



B. Green functions 

In addition to the lesser Green function defined in Eq. 
([5]), we define the following retarded Green functions: 

Gf^(a;,a;';i,t') = -*0(i-<')({V'a(2:,t),V'^(x',t')}) 

where the first one corresponds to degrees of freedom of 
the wire, while the second one corresponds to degrees of 
freedom of the wire and the j-th reservoir. 

The evaluation of these Green functions implies the 
solution of the Dyson equations. For the sake of sim- 
plicity in the notation, it is convenient to carry out the 
following gauge transformations -01^(2;) — V'/^(2;)e^''^*'^, 

xl- r iVj) ~^ e^^'^F^y^xl- r iVj)- Dyson equation for 
the retarded function reads 



{-idf ±kF-Hp± id,,}G^f,{x, x'- 1, t') - \bG^-^{x, x'; t, t')S{x' - xt) = 6{t - t')So.p + 

y^A, 2/°; - X,) +T. f dx"dt"G^^{x, x"- 1, x'; t", t') 

{-tdt,±k'^^ ±zdyJG^f,^{x,yf,t,t') = w,J2G^^{x,xf,t,my,-y^), 



(10) 
(11) 



where the upper and lower signs correspond to /3 — l,r 
and f3j — lj,rj, respectively, and I — r, r — I, while 
T,^^p{x" ,x';t" ,t') is the exact self-energy due to the in- 
teraction term with coupling constant g. 

Let us now notice that the operator 

-idt'±kF±idy>^ [gj^^iyj,y'^;t,t')} \ (12) 

is the inverse of the retarded Green function correspond- 
ing to the degrees of freedom j3j of the reservoir j. Thus, 



Eq. ([TT|) can be expressed as follows 

Gf^^.(x,2/°;t,i') = w.Y^I dt"G^f,{x,x,;t,t") 

xgf^iyly'^-XX). (13) 

Substituting the latter equation into Eq. pT|) and defin- 
ing 

^}'{x,x';t,t') = J2 S{x - Xj)S{x' - Xj)\wj\^ 

X5|.(y«,y^";t,t'), (14) 



4 



leads to 



{-tdf ±kF-tip± id,,)G%{x, x'; t, t') - XsG^-^ix, x'; t, t')5(x' - Xb) = 5{t - t')5c,p 

+ I dx''dt''G^^{x,x'';t,t''){J::;^{x'\x';t'\t') + ^''f^'{x'' ,x';t'\t')}. (15) 

I 

The lesser Green function entering the expression J dt"[A^{t,t")B<{t" ,t') + A< {t,t")B'^{t" ,t')], being the 

for the currents Ij can be obtained by means of advanced function [B'^]'^ = B^ . Thus 
Langreth rules from according to which given 
C«(t,t') = fdt"A"it,t")B^{t",t') then C<{t,t') = 



G<p^ y°; t') = w, / dt" G<Jx°, i, t")gj^ [y^ y°; t" , t') + G1(x°, x^ t, t")gl y°; t" , t') (16) 



where 17^ = [g^] ^ is the advanced Green function of the 
uncoupled reservoir. 

So far all the equations are exact. The crucial step to 
obtain the exact Green function by solving Dyson equa- 
tions is the evaluation of S™* , which corresponds to the 
fully dressed skeleton diagram for the self-energy corre- 
sponding to the electron-electron interaction, also taking 
into account the coupling to the two additional reservoirs 
as well as the backward impurity. We now introduce the 
following approximation for the limit of weak coupling 
to the reservoirs and the impurity such that Wj <C g and 
As <5: 



^7/3* (^i a;'; t') 



(17) 



where 



i:^^'\x,x';t,t') = [G«^^""(x,x';i,t')];;3' " 
{-idf ±kF±id^'}6^p, (18) 



is the self-energy of the infinite Luttinger wire with- 
out impurity and uncoupled from the reservoirs, while 
G^'^"**(a;, a;'; i') is the ensuing retarded Green func- 
tion. The approximation (|17p implies the evaluation of 
the self-energy associated to e-e interaction by neglecting 
vertex corrections due to the escape to the reservoirs and 
due to the scattering with the impurity. This approxi- 
mation is adequate only in the limit of small Wj and . 



Under this approximation in the e-e self-energy and 
performing a Fourier transform with respect to t — t' , Eq. 
(fT5| can be expressed as follows 



Gaf)ix, x']uj) = Gf^^"**(a;, x'; uj)+Y^ G^ai^. Xj,uj)Tr^^{xj,Xj\u:)Gf^''*\xj,x' ; uj)+G^Jx, Xb, cj)ABG^^"**(xf„ x'; uj) 

09) 



This equation allows for the evaluation of the retarded 
Green function. In what follows, we solve it at the lowest 
order in the backscattering term Xb and up to 0(w|), 
in the coupling to the voltage probes. We recall that 
ideal non-invasive probes correspond to keeping only up 
to 0{wj). It is important to notice that in the limit of 



vanishing Coulomb interaction (g — 0), the above equa- 
tion leads the exact retarded Green function of the prob- 
lem. 
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Currents flowing through the contact between the j'-th reservoir 

and the wire 



Substituting Eq. (jl6p in the definition of the current 
([7|), we get the foUowing exact equation for the current 



(20) 



Making use of the assumption of weak couphng be- 
tween the probes and the wire and weak amphtude 
in the back scattering term induced by the impurity, 
we evaluate the Green functions G^„(x°, a;°; w) and 



G^^{xj,x'j;uj) perturbatively up to 0{wj) and 0{Xb)- 
Concretely, this implies solving as 



j 

+G^„^""(x, Xb, ^)AbG^^""(x6, x'; lj), (21) 



while the lesser counterpart can be derived from ([21]) by where /j^'' oc /j^'' oc being 

recourse to Langreth rules (see above Eq. ([TB)) ). The 
explicit expression for G^^^"**(x, , w) is given in Ap- 
pendix A. After some algebra, the currents through the 
contacts can be expressed as follows 

/^.^/W-f/f , (22) 



iP = 



H E r ^\G<i''''\x,,x,,u:)g>{oo)-Gl;t-'\x,,x,,oj)g<{uj) 

+Ab[[G<„^""(2:, - X,;.:) (G^^""(a;, - x,;^))* + Gf^""(x, - a;)G#""(a:, - x,)]g>{.^) 
[G>„^""(a:, - x,;uj) (g#""(x, - x.-u:))* + Gf^""(x, - x,-u:)G^-'\x, - x,)]g<{u:)\ | 
^{Gf^«"(x„x„a.)G^^""(:.„a:„c.) [5>(c.).9<(a.) - g<(a.)g>(a.)] + 

G^/"" (a;j , , w)g> (w) - G^^^"" (a;, , a;^ , a;)g/ (cj) + 



(23) 



Aw] ^ 

i=l,2 



+ Gf^""(x„x„a.)gf(a.) 



(24) 



+ Gi'^^-'\x.,x„u:)gt{u^) [g</^""(x„ x., ^)5>(c.) - G>/""(x„ a;,, c.).9<(^)] } 



The term /j^^ corresponds to the limit of ideal non- by dropping the second term in Eq. (|2T|) and the en- 
invasive contacts considered in Ref. [l^. It is derived 
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suing term in the lesser counterpart. This leads to the 
exact solution 0{wj) of ^5)) and for — 0. In 
the second-order solutions ([21]) we have introduced the 
additional approximation of neglecting vertex corrections 
oc Wj and oc in the evaluation of the many-body self- 
energy S*"*. Notice that the two probes are completely 
uncorrelated within the "non-invasive" component 
In the higher order contribution I^^^^ it is possible to dis- 
tinguish two kinds of terms. On one hand, those oc w*, 
account for the effect of dephasing and resistive behavior 
induced by the inelastic scattering processes due to the 
coupling to the reservoirs. On the other hand, terms 
oc wlw2 describe interference effects between the two 
probes. 

It is now convenient to express the lesser and greater 
Green function in terms of spectral functions: 

g,<,>,L«tt(-^^ = A<'>(w)p„(a;i -a;2;w), 

g<'>ic.) = A<->Hp,(c.) (25) 



with A< ., (w) = inF,aj(w), A> ,,-(w) = - nF,aj(w)], 
being np.a 

tion where the upper and lower signs corresponds, re- 
spectively to the right and left movers of the wire, and 
nFj(w) = l/(e("~'^j)/^ + 1), fj,j being the chemical 
potentials of the electrons in the j reservoir, relative 
to the mean chemical potential /i of the wire. T is 
the temperature, which we assume to be the same for 
the wire and the probes, while Pa{i^) = «G'^^^""(a;i — 
X2;oj) — i[Gf^^""(a;2 — xi;uj)]* is the spectral function 
for the a movers in the Luttinger model, and Pj{u)) = 
—21m[g^{(jj)\, is the spectral density of the j probe. Re- 
placing in (|24|) and ([25]) . the full expression for the cur- 
rent reads: 



duj 
2^ 



np^ait^) - nF,j{u!)\pj{uj){pa[0,uj) 1 -I- 2w|Re {g^'^"^*^ [G,uj)gf{Lo) 



-t-2ABRe 



-2«;>| 



E 



2w|Re 



Pa{xj - x--uj)*Gf^'^*\xj - xy,uj)gf{uj) } 



duj 
2^ 



where we use the notation j,j such that 1 = 2 and 2 = 1. 



D. Voltage drop and four-terminal resistance 

The chemical potentials pj in the expressions of the 
previous subsection must be set to satisfy the condition of 
local electrochemical equilibrium between the probes and 
the wire. This implies vanishing flows Ij = 0, j — 1,2, 
with the currents defined in Eq. , and the two chem- 
ical potentials pj must satisfy these constraints. In the 
case of non-invasive probes, the two probes are com- 
pletely uncorrelated, and the problem can be reduced to 
that of the wire coupled to a single probe, which senses 
the local chemical potential of the wire. Instead, in the 
present case, we have to solve a system of two non-linear 
equations to calculate pi and p2, from where the volt- 
age drop AV^ — Pi — P2 between the points xi and X2 
of the wire coupled to the two probes can be evaluated. 
This voltage drop contains not only information of the 
scattering processes in the wire that are independent of 
the coupling to the probes, but also of inelastic scatter- 
ing processes and interference effects introduced by the 
probes themselves. The four-terminal resistance can be 



evaluated from Eq. ([T]) and the ratio between the four- 
terminal and two terminal resistance i?2t — V/I results 



R4t 
R2t 



AV 



(27) 



The two chemical potentials are evaluated numerically 
from Eq. (j26p . with the Green functions given in Ap- 
pendix A. 



III. RESULTS WITHOUT IMPURITY 

In this section we show results for the ratio between 
Rit/R2t in the case of A^ = 0. It is important to men- 
tion that in the limit of non-invasive probes, this ratio 
vanishes identically under this case, and all the features 
in the behavior of the resistance discussed in this section 
are solely due to the invasive nature of the probes. 

We characterize the strength of e-e interactions with 
the parameter K — (1 + ^)^^/^. The limit of non- 
interacting electrons corresponds Xo K — 1 while typi- 
cal values of K (experimentally determined in transport 
measures in nanodevices^^) are in the range 0.25 < K < 
0.75. 



FIG. 2: (Color online) i?4t/_R2t as a function of the voltage 
V, for different values of the e-e interaction strength K. The 
positions of the probes are xi = —10 and X2 = 10, and the 
strength of the couplings are wi — W2 ~ 0.1. The tempera- 
ture is T = 0. 




FIG. 3: (Color online) RAt/R2t as a function of the position 
of the second probe X2, given the first probe fixed at x\ — 
— 10.. The strength of the e-e interaction is K — 0.7, and the 
couplings ^11=102 = 0.1. 



FIG. 5: (Color online) _R4t/i?2t as a function of the voltage V 
and the temperature T, for K — 0.7. The positions of the 
probes are xi = —10. end X2 ~ 10., and the strength of the 
couplings are wi = W2 = 0.1. 




FIG. 6: (Color online) ii4t/-R2t as a function of the temper- 
ature T, for K = 0.7 and three values of V. The positions 
of the probes are xi = —10. end X2 — 10., and the strength 
of the couplings are uii = W2 = 0.1. Inset: low temperature 
regime. 
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FIG. 4: (Color online) Ru/R2t as a function of the position of 
the second probe X2, given the first probe fixed at xi = —10., 
for different values of the e-e interaction strength K. We set 
T = Q,wi^W2^ 0.1 and V = 1. 



Results for R4t/R2t as a function of the bias voltage 
V, for the probes connected at two fixed positions and 
different values of the e-e interaction K are shown in Fig. 
[5] In order to gain insight on the behavior of the ratio 
between resistances, let us notice that for vanishing bias 
voltage V, the voltage drop AV^ and thus R4t/R2t should 
be also vanishing. It is, therefore, not surprising that for 
low enough V, Ru/ R2t displays a power law behavior as 
a function of V, 

^ oc y27+i. (28) 

R2t 

The exponent 7 is related to the Luttinger parameter as 
= (K + - 2)/4. The latter result can be rather 
straightforwardly derived from an expansion of Ij for low 
V. On the other hand, a classical ohmic-like resistive be- 
havior implies a constant value of R4t/R2t- In Fig. [21 
it can be seen that such a behavior is approximately at- 
tained when the bias voltage overcomes a value Vc, which 
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satisfies eVc « hv/{2K | zi — 12 This energy scale 
can be understood by noticing that v/K is the plasmon 
velocity along the wire and Tp — 2K \ xi — X2 \ /v is 
the time that these excitations take for a round trip be- 
tween the probes. The latter defines the characteristic 
time for the inelastic back-scattering processes. Notice, 
that although the Luttinger wire is an elastic system, 
where electrons propagate ballistically, the coupled volt- 
age probes act as a dissipative bath. In fact, it is precisely 
the coupling to reservoirs the mechanism usuall y fo llowed 
in the literature (see, for example Refs. [13 andl25[) in or- 
der to model Ohmic dissipation. In the present case, 
assuming that the bias is applied from left to right, the 
Fermi energy of right-moving electrons is an amount eV 
higher than that of the left-moving ones. Then, the en- 
ergy associated to the crossover voltage eVc corresponds 
to the energy dissipated in the contacts, in a process in 
which an electron with the Fermi energy travels with 
velocity v/K from the left probe, connected at xi, to the 
right one, connected at X2, it is backscattered at X2 and 
comes back to xi with a Fermi energy fj.^. An estimate 
for the energy transfer involved in the dissipative process 
is, precisely, Tit~^. The above argument can be easily 
reconstructed for the case of a bias with opposite sign, 
in which case, the energy is inelastically transferred from 
right to left movers. Notice that in any case, the energy 
dissipated at the contacts is associated to a voltage drop 
that has the same direction as the external bias, as is ex- 
pected for a classical Ohmic-like process. Interestingly, 
2t~^ is equivalent to the ballistic frequency defined in 
Ref. [19I for an interacting Luttinger wire of finite length 
connected to reservoirs. 

To summarize, for a given separation \ Xi ~ X2 \ be- 
tween the two probes, Vc defines the crossover voltage 
for which inelastic back-scattering processes between the 
two points become active. Notice that the low voltage 
regime F < Vc so defined, depends on the e-e interac- 
tion strength 5, being wider for stronger g (lower K). In 
general, the effect of this interaction is to decrease the 
resistance. A closer analysis of Fig. [2] for V > Vc reveals 
that i?4t /i?2t as a function of the bias voltage V displays 
oscillatory features. This can be naturally interpreted as 
the consequence of interference effects between the two 
probes. From fits of the numerical data, we found that 
they can be very well reproduced by a function of the 
form: 

^^A + B sin{KV{x2 ~ xi))/V^^+\ V > Vc, (29) 

with A and B depending on K while proportional to 
w"^ , although we have not derived this result analytically 
from Eq. (j26p . It is anyway interesting that a similar 
resistive behavior is obtained in a Luttinger wire of finite 
length in the presence of back-scattering processes (see 
Refs. [19112^ . Another interesting observation is that 
the saturation value A decreases for increasing electron- 
electron interactions, This indicates that the latter tend 
to screen the inelastic scattering processes introduced by 
the coupling to the probes. 



In Fig. |3]we show the behavior of the ratio between re- 
sistances with the position of one of the probes kept fixed 
while the position of the second one is moved along the 
wire. This pattern reveals that the functional behavior 
is 

^ ^ A + Bsm{2kF{x2-xi))sm{KV{x2^xi))/V^'<+\ 
n2t 

(30) 

within the high V > Vc voltage regime, corresponding, 
respectively, to solid and dashed lines in the Figure. The 
2kF modulation resembles the behavior found in the volt- 
age profile of non-invasive probes in a system with an im- 
purity, which is observed both in non-interactingiSiiiii^ 
and interacting systema^S. In those cases the origin 
is the occurrence of interference in the electronic wave 
packet generated by the back-scattering processes that 
take place at the impurity. In the present case, the 
interference is originated by scattering processes at the 
probes. Unlike the behavior for non-invasive probes, in 
our case the voltage drop induced by the probes has the 
same sign as the applied external voltage. This means 
that the four-terminal resistance for invasive probes in 
a clean wire is always positive, in spite of the Friedel- 
like 2k F oscillations. This is in strong contrast to the 
case of non-invasive probes, where these oscillations pro- 
vide a mechanism for i?4t < 0. Fig. 2] illustrates the 
same situation but for fixed voltage and varying K . One 
sees that, in general, larger values of the e-e interactions 
produce smaller values of i?4t/i?2t- Then we conclude 
that, although one cannot have negative values of the 
four terminal resistance in the absence of impurities, e-e 
interactions tend to facilitate that possibility. 

In Fig. [5] we show the effect of the temperature in 
the behavior of i?4t/i?2t- It is clear that, as the tem- 
perature increases, the oscillations discussed in Fig. [2] 
within the high voltage regime, tend to be wiped out and 
the resistance evolves to a constant value. This behavior 
is depicted in more detail in Fig. [SJ where we display 
Rit /R2t as function of T for three different values of the 
bias voltage V. In analogy with the previously discussed 
behavior found at T = 0, as function of V (Figl2]), there 
is a crossover temperature Tc w hv/(2K \ X1—X2 |) which 
allows to distinguish low and high temperature regimes. 
For low temperatures (T < Tc) , we have verified that the 
ratio between resistances behaves as 

^c^a + bT^^+\ (31) 

where a and b depend on V. For high temperatures 
RAt/R2t tends to a constant value. As the temperature 
increases, coherence tends to disappear. For this reason, 
no signature of the oscillatory behavior observed in Fig. 
[5] is found here. The interplay between T and V gives 
rise to the possible occurrence of a global maximum of 
the i?4t / i?2t • The value of temperature that corresponds 
to this maximum, when it is present, depends on K and 

\xi- X2 |. 



9 




-15 -10 -5 5 10 15 5 10 15 20 25 30 35 40 



FIG. 7: (Color online) Ru/Rit as a function of the position 
of the impurity Xfj, for probes connected at x\ = —10. and 
X2 ~ 10. The bias voltage is V = 0.5, the strength of the 
e-e interaction is K = 0.7, and the couplings wi — W2 ~ 
0.1. Different plots correspond to different strengths of the 
backscattering term representing the impurity. 



IV. RESULTS WITH IMPURITY 

In this section we analyze the behavior of Ru/R2t at 
T — 0, for a wire with an impurity with backscattering 
strength A^. In the case of non-invasive probes, the lo- 
cal voltage displays 2kF Friedel-like oscillations with con- 
stant amplitude for non-interacting electron o^°i^^'^'^ , and 
with modulated amplitude in the case of and interacting 
wirei^ 

Figure [7] shows R/^t / R2t for the probes connected at 
fixed positions, as a function of the position of the im- 
purity Xb- Friedel-like oscillations with period 2kp are 
identified, with an increasing amplitude for increasing 
back-scattering strength. As in the case of non-invasive 
probes, the amplitude is modulated for interacting elec- 
trons, the local voltage achieving the highest amplitudes 
at the position of the impurity. Unlike the case of non- 
invasive probes, the oscillations take place with respect to 
a constant non-vanishing value, which is determined by 
the degree of coupling of the probes. For the parameters 
shown in the figure, Ra is always a positive quantity. 

Besides interference effects, it is clear that the cou- 
pling of the probes generates classical resistive behavior 
through inelastic scattering processes, while the elastic 
scattering induced by the impurity induce Friedel oscil- 
lations. The first type of processes takes place with a 
strength oc w'j /A', where A' is the bandwidth of the reser- 
voirs, while the second one takes place with a strength 
Xb- The two mechanisms are competitive regarding the 
possibility of having i?4t < 0. In Fig. |8] we analyze, 
precisely, this possibility. To this end, we have fixed 
the first probe at the position xi = xi, = —10, where 
the minimum R4t/R2t is achieved, considering different 
positions for the second probe X2- For each of these 
configurations we then vary the ratio Wj/Xb, to define 
{wj / XB)iim, at which Rit/R2t = 0. The corresponding 



FIG. 8: (Color online) The ratio {w^ /XB)iim at which Ra = 
for x\ = xi, = —10, as a function of the position of the 
second probe X2- Different plots correspond to different e-e 
interaction parameters K. The applied voltage is 1/ = 0.5. 

results are plotted in the figure for different e-e inter- 
actions. For w'^/Xb > {wj /XB)iim the ratio Rit/R2t is 
positive for any value of Xh. On the other hand, the condi- 
tion Wj/Xb < {w-j /XB)hm defines the values of coupling 
strength for which a negative four terminal resistance is 
possible, depending on the position of the impurity. 

A very interesting and subtle issue that is also revealed 
by our analysis concerns the role of e-e interactions in 
the possible occurrence of a negative four-terminal re- 
sistance. Based on the results obtained for non-invasive 
probesiSiiiii^iSS one would expect that e-e interactions 
oppose to such possibility, owing to the fact that for 
stronger interactions (smaller values of K) the ampli- 
tude of the oscillations coming from the presence of the 
impurity diminishes. However, in the present case this 
effect competes with the global "upward" shift coming 
from the contribution of Jj . In other words, as already 
pointed out in Section III, the weak invasiveness of the 
probes, which in our formulation is contained in /j pro- 
duces a voltage drop that has the same sign of the bias 
V. It turns out that the magnitude of such a shift also 
depends on K, and it decreases for increasing interac- 
tions (decreasing K), as shown in Figures ^ and 
The combination of these two effects gives rise to the 
result depicted in Fig. ([81), where one sees that for suf- 
ficiently separated probes, e-e interactions facilitate the 
occurrence of a negative four-terminal resistance. 

V. SUMMARY AND CONCLUSIONS 

We have analyzed the behavior of the four terminal re- 
sistance in a biased quantum wire with an impurity. We 
have modeled the wire by an infinite-length Luttinger 
wire where the bias voltage is represented by different 
chemical potentials for the left and right movers, and the 
impurity by a backscattering term. We have also intro- 
duced models for the probes, which consist in reservoirs 
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of non-interacting electrons. These systems are locally 
weakly coupled to the wire and have chemical potentials 
satisfying the conditions of vanishing electronic currents 
between the reservoirs and the wire. The difference be- 
tween the so determined chemical potentials defines the 
voltage drop, from where the ratio between the four- 
terminal and two-terminal resistance can be calculated. 
We have solved the problem within perturbation theory 
in the impurity strength and the tunneling parameter 
defining the coupling between the probes and the wire, 
within the framework of non-equilibrium Green functions 
formalism. We have neglected vertex corrections in the 
self-energy for the e-e interaction associated to inelas- 
tic scattering processes due to the escape to the leads 
and elastic scattering processes at the impurity. Since 
we have assumed that these two parameters are small 
enough, the latter is expected to be a reliable approxi- 
mation. 

We have analyzed the voltage drop beyond the non- 
invasive assumption for the coupling of the probes to 
the wire. That is, we have studied, not only the volt- 
age drop originated by elastic scattering processes along 
the wire, but also the effects introduced by the coupling 
to the probes, itself. We have shown that the inelastic 
scattering processes due to the invasive coupling of the 
probes induce a voltage drop with a power law behavior 
as a function of the bias voltage for low values of this 
parameter, with an exponent determined by the e-e in- 
teraction. In the limit of non-interacting electrons, this 
reduces to a linear dependence as a function of the bias 
voltage. This behavior has classical and quantum fea- 
tures, since the voltage drop is always in the same sense 
of the applied voltage but displays a pattern of oscilla- 
tions indicating quantum interference between the two 
probes. These features, are, however, screened as the e-e 
interaction increases. In our calculations, we have consid- 
ered an infinite wire. However, the separation between 
the probes sets a natural length scale in the problem, 
which determines the crossover value of the bias voltage 
for which inelastic scattering processes become active. In 
the case of an interacting wire with finite length, we ex- 
pect that our results remain valid provided the length 
of the wire is much larger than the separation between 
the probes. In the presence of an impurity, the elastic 
backward scattering processes and 2kp oscillations de- 
tected by non-invasive probes^ are superimposed to the 
inelastic processes introduced by the probes. 

Our results have an important outcome in relation to 
experimental measurements of four-terminal resistance in 
real systems. That is, for invasive probes, elastic effects 
like those generated by backscattering processes by impu- 
rities can still lead to a voltage drop that opposes to the 
applied voltage, giving rise to a negative four-terminal 
resistance. However, the amplitude for these processes 
must be high enough in order to overcome the classical 
resistive effect introduced by the probes. 

As far as the e-e interaction effects are concerned, 
they play a fundamental role in the calculated magni- 
tudes. For higher e-e interaction, the oscillations ampli- 
tude coming from the impurity decreases. The ampli- 



tude of the global shift coming from the interaction of 
the probes also decreases for stronger interactions. We 
have shown that if the separation of the probes is large 
enough, the possibility of measuring a negative resistance 
increases for stronger interactions. 



VI. ACKNOWLEDGMENTS 

We acknowledge support from CONICET and AN- 
PCYT, Argentina, as well as UBACYT and the J. S. 
Guggenheim Memorial Foundation (LA). 



Appendix A: Green functions and spectral functions 
for the Luttinger wire and the reservoirs 

We can follow the procedure of Ref. [3l| to evalu- 
ate the retarded Green functions of the Luttinger wire 
and calculate the spectral density from pa{x,x';u!) = 
■^QR.Lutt^^ _ _ QR.,Lutt(^^i _ x-uj)*]. The rcsuh 

is 



Pa{x,x' ]UJ + Pa) — C exp 



kp] X 

V 



1 27 



X(/)(7,27+ l,±2ixLj/'u), (Al) 

where (j) is Kummer's Hypergeometric function, 7 = {K+ 
— 2)/4 and v — 1/K. In order to perform numerical 
calculations, we introduce an energy cutoff A such that 
p{x, UJ + p) ^ p{x,uj + p)&{A — \uj\). C is a normalization 
constant, which can be evaluated by the normalization 
condition 



duj I 



= 1 



(A2) 



The retarded and advanced green function are then cal- 
culated using the Kramers-Kronig relations 



G 



R.Lutt 



{x,uj) 



did' p{x,uj') 
27r UJ ^ uj' + irj 

duj' Ile[p{x, Lu')] + ilm[p{x, u') 



2tt 



uj — lu' + irj 



{A3) 



f duj' p{x,uj') i , , , . ^, 

The real part is evaluated numerically by using the pro- 
cedure explained in Ref. [s^. We have verified that with 
a cutoff A w 20 the evaluated voltage drop AV^ is inde- 
pendent of this cutoff. 

For the reservoirs, we consider a constant density of 
states within a cutoff ±A'. So, the retarded Green func- 
tion for the probes can be calculated using the Kramers- 
Kronig relations and gives 



gf{Q,uj + p) 



2A 



1 , i^ + A', 
-ln|- 



'w- A' 



2A 



-e(A' 



(A5) 
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